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EXISTENCE  AND  UNIQUENESS  THEOREMS  FOR  THE 
NEUTRON  TRANSPORT  EQUATION 

Abstract 

In  an  attempt  to  understand  the  conditions  under  which  the  neutron  trans¬ 
port  equation  lias  solutions,  and  the  properties  of  those  solutions,  a  number  of 
existence  and  uniqueness  theorems  are  proved.  One  finds  that  the  properties  of 
the  solution  are  closely  related  to  the  boundedness  of  the  source  as  well  as  to 
certain  velocity-space  Integrals  of  the  scattering  kernel.  Both  time -dependent 
and  time-independent  equations  are  considered  as  are  also  the  time -dependent  and 
time -Independent  adjoint  equations.  Although  only  a  very  few  of  all  possible 
existence  and  uniqueness  theorems  for  these  equations  are  considered  here,  the 
work  may  serve  as  a  guide  to  the  treatment  of  similar  problems. 
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I.  Introduction 


Various  theorems  concerning  the  existence  and  uniqueness  of  solution  to 
the  neutron  transport  equation  have  appeared  in  the  literature.  For  example. 

Case  (l)  has  proved  uniqueness  for  the  one-speed,  time-dependent  equation 
under  the  assumption  that  the  kernel  is  rotationally  invariant.  He  has,  in 
addition,  shown  that  under  the  same  assumptions  the  solution  of  the  time- 
independent  equation  is  imique  if  o(x)>  'the  mean  number  of  neutrons  emitted 
per  collision,  is  everj'where  less  than  one.  Olhoeft  (2)  has  considered  the 
more  general  velocity-dependent  case,  and  has  shown,  subject  to  the  same 
restriction  on  c,  that  a  unique,  integrable  solution  exists  for  the  time- 
independent  equation.  Davison  (3)  has  made  some  rather  general  remarks  con¬ 
cerning  existence  and  uniqueness  for  the  time -dependent  case,  but  has  only 
outlined  the  methods  of  proof  and  has  actually  said  little  about  the  restrictions 
which  must  be  Imposed  in  order  that  the  theorems  be  true. 

Basically,  the  situation  is  the  following:  There  are  a  number  of  pos¬ 
sible  restrictions  which  one  can  imagine  might  be  applied  to  the  cross-sections 
and  sources  appearing  in  the  neutron  treuisport  equation.  For  certain  of  these 
restrictions  it  is  possible  to  prove  that  continuous  solutions  exist;  for  other 
sets  of  restrictions  the  solutions  may  not  be  continuous  but  still  integrable 
functions  of  one  or  more  of  the  Independent  variables  involved  (i.e.,  position, 
velocity,  and  [in  the  time -dependent  case]  time).  For  other  restrictions,  it 
may  be  Impossible  to  prove  anything.  We  have  investigated  a  large  number  of 
possible  restrictions  which  might  reasonably  be  imposed  upon  the  cross-sections 
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and  sources,  and  investigated  the  existence  of  unique  solutions  for  each  case. 

In  this  way,  we  have  tried  to  bring  some  order  Into  the  chaos  of  "obvious" 
or  partially  proved  results  which  at  the  present  exist  in  the  literature. 

In  addition,  we  consider  not  only  the  transport  equation  but  the  time- 
dependent  and  time-independent  "adjoint  equations,"  and  investigate  the  suf¬ 
ficient  conditions  that  unique  solutions  exist  for  those  equations.  It  turns 
out  that  there  are  many  cases  in  which  existence  and  uniqueness  theorems  can 
be  proved  for  one  or  the  other,  but  not  both. 

In  Sec.  II  of  this  paper,  we  convert  the  transport  and  adjoint  equations  to 
integral  equations  in  the  usual  manner.  Then,  in  Sec.  Ill,  we  discuss  the  re¬ 
strictions  on  the  sources  and  cross-sections  which  will,  for  physical  reasons, 
be  applied  In  ail  cases. 

Then,  In  Secs.  IV  and  V  we  consider  various  existence  and  uniqueness  theorems 
for  the  time-dependent  transport  and  adjoint  equations.  In  Sec.  VI  we  consider 
theorems  for  the  time -independent  equations  (both  transport  and  adjoint)  and 
finally,  in  Sec.  VII,  we  discuss  certain  "by-products"  of  the  theorems  --  a 
formula  for  the  minimum  critical  size  of  a  reactor  and  a  proof  that  the  time- 
dependence  of  the  solutions  of  the  time -dependent  equations  must  obey  certain 
restrictions. 

We  have  made  no  attempt  to  be  comprehensive  in  our  treatment  of  the  existence 
and  uniqueness  problem.  Indeed,  one  cein  imagine  virtually  an  infinite  number 
of  possible  theorems  similar  to  those  which  we  prove.  However,  we  have  pre¬ 
sented  some  of  those  which  we  feel  have  the  most  intrinsic  interest  and,  in  ad¬ 
dition,  demonstrate  well  the  method  of  proof.  In  this  way,  if  any  of  the 
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conditions  vhich  we  consider  are  not  met  in  a  particular  problem,  our  work  may 
serve  as  a  guide  to  the  correct  treatment. 

II.  The  Integral  Equation  Formulation 


A.  The  Transport  Equations 

The  time -dependent  transport  equation  is 


di|r(r,v,t) 

ht 


+  v(n*7+a(£>v))\|f(r,v,t) 


q(r,v,t) 


+ 


'a(r,y '■*v)\|r(r,v' ,t)d3v  , 


(la) 


while,  in  the  stationary  limit,  we  have 


v(n*y+cT(r,y))\|f(r,y) 


cr(r,y'-*y)i|r(r,y '  )d®v'  . 


(lb) 


Here  \|r  is  the  neutron  ang\ilar  density  (l.e.,  the  one-particle  distribution  function); 
y  =  vn  is  the  neutron  velocity;  a(r,y)  is  the  total  cross  section,  and  the  kernel 
a(r,y'-*y)  is  the  cross  section  for  a  neutron  of  velocity  y'  to  be  emitted  into  d^ 
about  y. 

Eq.  (l)  may  be  converted  into  integral  equations  in  the  usvial  way,  i.e.,  by 
introducing  the  Green's  function  of  the  left  hand  side.  This  Green's  function  is 
well  known  (4).  We  obtain  for  the  Integral  eqiiatlon  equivalent  to  (la) 

\|r(r,y,t)  =  Q(r,y,t)  +  ^  dt 

X  \jf(r-y(t-t''),v’,t')exp 


va(r-v(t-t"),v)dt" 


(2a) 


'J  d3v'v'cj(r-y(t-t' ),  y'^^y) 
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where 


Q(r,v,t) 


ilf(r-vt,.  V,  0)exp 


va( r-v( t-t ' ) , v)dt ' 


t-t' ) ,v,t' )exp 


t 

va(r-v(t-t"),v)dt" 


(2b) 


Actually,  we  shall  prove  theorems  involving  the  existence  and  uniqueness  of 
solutions  within  a  given  volume  of  space  V  boxmded  by  a  surface  S  when  the  incoming 
(or,  in  the  case  of  the  adjoint  equations,  the  outgoing)  angular  distribution  is 
specified  on  S.  The  simplest  way  to  modify  Eq.  (2)  to  describe  this  situation  is 
to  define  q(r,v,t),  \|r(r,v,0),  a(r,v) ,  and  a(r,v'+v)  to  vanish  for  r  ^  V.  Then 
the  specified  Incident  angular  dlstrlbutioni|lj:(rg,v,t)  is  replaced  by  a  surface 
source  <ls(fs>Z>^)  ®  usxial  manner  (k). 


■s  0  otherwise. 


(5  ) 


where  Hq  is  the  outward  normal  to  S.  Then  Eq.  (2)  still  applies  (understanding 
the  redefinitions  of  q,  \|f,  a,  and  c(r,v'-*;v)  mentioned  above)  with  an  additional 
term  Qs  added  to  Q,  Eq.  (2b)  to  account  for  the  surface  source: 


i  -0((r,r-Rsn,v) 


(4  ) 


Here  Of  is  the  usual  optical  path 


a(r,v',v)  =  ^ 


d8a(r-a=. )  , 

■r 


(5  ) 


where  Rg  is  the  distance  frcan  r  to  the  surface  S  along  the  direction  -fl  (See  Fig. 
l).  Note!  Rg  ■  Rg(r^n) . 


Eq.  (2),  thus  modified,  is  the  general  integral  equation  formulation  of 
the  time  dependent  neutron  transport  problem  with  which  we  shall  be  concerned. 
It  is  convenient  to  rewrite  it  in  the  somewhat  more  tractable  form 


t(r,v,t)  =  Q'(r,v,t)  +  Kt(r,v,t), 


(6a) 


where  K  is  the  integral  operator 

•  t 


Kf(r,v,t)  =  ^  dt'J  6.^r'J^i^v'b(T'-T+v(t-t'))exp  -J  va(r-v(t-t")  ,v)dt" 


X  v'a(r ',v'-^)f(\r ' ,v' ,t' )  , 


(6b) 


and' 


Q'(r,v,t)  -  Q(r,v,t)  +  Qs(r,v,t) 


(6c) 


The  Integral  equation  fonnulatlon  of  the  time-independent  transport  equation 
(lb)  is  obtained  in  just  the  same  way,  i.e.,  by  introducing  the  Green's  function 
of  the  left  side  of  Eq.  (lb).  We  obtain 
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where  A  Is  the  integral  operator 


Af(r,v)  *  ^  dR^ d%'^d®v'6(r'-r+Rn)exp  fa(r,r' ,v)| 


L 


v'a(r'  ,v'-»v)f(r'  .v' ),  (7b) 


Q'(r,v)  =  <ls(£s»v)exp  ^0((r,r-Rgn,v^  q(r-Rft,v)exp  j^a(r,r-Rfl,v^  dR,  (To) 


and 


♦(rjv)  =  vt(r,v)  .  (7d) 

Again  we  mention  that  the  cross-sections  and  sources  have  been  defined  to  vanish 
for  r  1^.  V  (or  for  t  <  O). 

B.  The  Adjoint  Equations 

The  time  dependent  adjoint  equation  is  defined  to  be 


^(r,v>t) 

9t 


+  v(-n' V+o(r,v))t(r,v,t)  =  q(r,v,t) 


vj' d\fa(r,r^' )\|f(r,v' ) , 


(8a) 


while  the  time-independent  adjoint  equation  is 


The  adjoint  of  a  solution  of  Eqs.  (la,  lb)  will  be  taken  to  be  the  solution  of 
Eqs.  (8a,  8b)  subject  to  appropriate  boundary  conditions. 

These  boundary  conditions  are  that  the  outgoing  density  on  S,  to(£s>Z>'*')  > 
rather  than  the  incoming,  will  be  specified.  As  before  we  replace  the  boundary 
conditions  by  a  surface  source 

=  v|n*nQ|\|fr(ls>Z>'t)>  n*no>0,  t  >  0, 

=  0  otherwise  .  (  9  ) 


We  now  convert  Eq  (8)  into  integral  eqtaations  just  as  for  the  case  of  trans¬ 
port  equations, obtaining  for  the  time -dependent  equation 

iif(r,v,t)  =  Q'(r,v,t)  +  Ki|/(r,v  t),  (lOa) 

where  K  is  the  integral  operator 


Kf(r,v,t)  =  dt'y  d\*y'd\'B(r'-r-v(t-t' ) )exp 


va(r+v(t-t") ,y)dt' 


X  a(r' ,v^' )f  (r' ,v' ,t' ), 


(lOb) 


and 


Q'(r,v,t)  =  \|f(r+vt,v,o)exp 


va(r+v(t-t") ,v)dt 


r.t 

/  dt'q(r+y(t-t' ),y,t)exp 

-/  va(r+v(t-t"),v)dt’ 
^t’ 

+  i  qs(rs>v,t-Rs/v)exp 

V  “  “ 


-0((r,r+R8i?ij) 


(10c) 
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As  tefore  cross-sections  and  the  initial  distribution  are  defined  to  vanish 


for  r  4  ¥•  Also,  Rg  «  Rs(£>i})  distance  from  r  to  the  surface 
direction  +n  (rather  than  along  the  direction  -n  as  in  the  case  of  the 
equation).  See  Fig.  2. 


For  the  time-independent  adjoint  equation  we  obtain  similarly 

=  -I  5'(r,v)  +  A^((r,v)  , 


where  A  is  defined  by 


Af(r,v) 


dl^ d\'6(r'-r-Rn)exp 

a(r',:^')f(r',v'). 


-a(r,r',v) 


and 


Q'(r,v)  =  qg(rg,v)exp 


'-a(r£-i^,v) 


r\ 

- 

/  q(r+Rfi,v)exp 

-a(r,r+Rfi,v) 

dR. 


along  the 
transport 


(11a) 


(11b) 


(11c) 


8 


III.  Basic  Definitions  and  Assiamptlons 


In  all  of  our  suteequent  discussions  we  shall  assume  that  the  source  functions 
Q'(r,v,t)  and  Q'(r,v,t)  as  well  as  the  cross-sections  and  a(r,v)  obey 

certain  physically  reasonable  conditions. 

1.  Q'  and  5'  are  positive  for  all  values  of  their  arguments. 

2.  Either  Q'  and  Q'  are  bounded  or  they  can  be  written  as  the  product 
of  bounded  functions  multiplied  by  delta  functions  in  one  or  more  of 
their  argimients. 

3.  There  exists  a  Vq  <  “  such  that  for  |v|  >  Vq,  Q'  and  Q'  vanish  identi¬ 
cally. 

4.  There  exists  a  vi<  »  such  that  a(r,vW)  =  0,  for  v  >v'  >  va,. 

Assumptions  (3)  and  (U)  permit  us  to  avoid  any  difficulties  that  the  Infinite 

range  of  the  velocity  variable  might  otherwise  introduce  since  together  they  imply 
that  there  are  no  neutrons  present  with  speeds  greater  than  Vjj  =  Max(vQ,vi).  Then 
the  Integrals  over  v'  in  any  of  the  Integral  equations  can  be  written  as 


I 


d^v'  -*• 


dn 


vm 

V'2dv' 


(12) 


These  restrictions  can  actually  be  relaxed  to  some  extent;  we  might  assume  that 
Q'  and  Q'  go  to  zero  sufficiently  rapidly  as  v  -*•  »  that,  if  other  restrictions 
are  obeyed,  the  infinite  range  of  v  need  not  cause  any  trouble.  However,  assump¬ 
tions  (3)  and  (4)  are  physically  reasonable,  and  so  we  shall  make  no  attempts  to 
relax  them. 
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5*  o(r,v'-^)  can  be  written  in  the  form 


a(r,v'-^)  »  ^Ni(r)oi{v'-*v)  , 


(13) 


where  Ni(r)  is  bounded.  (Actually  %(r)  represents  a  density  of  nuclei, 
and  oi  is  a  microscopic  cross-section,  so  we  are  merely  assuming  that 
there  are  no  infinite  concentrations  of  atoms  present  in  the  systems 
that  we  consider.) 

6.  The  cross-section  a(r,v)  can  be  written  in  the  same  form: 


r(r,v)  = 


^Ni(r)ai(v)  . 


(U) 


7.  vai(v)  is  boionded.  We  expect  oi^(v)  to  be  bounded  except  possibly 
for  v-K),  in  which  case  we  admit  ai{v)~  l/v. 

8.  oi(r,v),  c(r,v'-*2),  ai(v)  and  (Ji(v'-*v)  are  all  positive. 

9.  The  function 

c(r,v)  =  /q(£>rz')d3v' 
a(r,v) 

is  positive.  This  assumption  is  also  physically  reasonable,  since 
c(r,v)  represents  the  mean  number  of  neutrons  emitted  per  collision. 

We  define  a  similar  microscopic  quantity  Ij^(y)  which  is  also  bounded 

jCi(v^')d®v' 


(15) 


£i(l) 


ai(v) 


(16) 
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We  shall  also  have  occasion  to  use  two  further  functions, 


c'(r,v) 


•*v)d®v' 


a(r,v) 


and 


c"(r,v) 


o(r,y'->y)d^v* 

a(r,y) 


(iTa) 


(17b) 


It  may  be  noted  that  o'  and  c"  may  not  always  be  bounded  but  they  (like  c) 
are  always  positive.  Similarly  t'  a  functions 


and 


5i'H 


5l"(v) 


ai(y '-►y)d®v' 
Ol(l) 


ai(y'-*y)d3v' 

ai(y) 


may  not  be  bounded  but,  like  li(y),  they  are  positive. 


(I8e) 


(iSb) 


IV.  The  Time-Dependent  Transport  Equation  . 

Bearing  In  mind  the  restrictions  discussed  In  Sec.  Ill,  we  consider  now 
various  existence  and  uniqueness  theorems  for  the  time -dependent  transport  equation. 
Theorem  1.  Let  Q'(r,y,t)  be  bounded.  Then  If  c"(r,y)  Is  bounded,  a  unique, 
positive,  and  continuous  solution. .to  the  time -dependent  transport  equation 
exists. 
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W«  prove  the  theorem  by  constructing  the  Neuram  series  solution  to  Eq.  (6)t 


tCrjVjt) 


00 

n=o 


where 


\|ro(r,v,t)  =  Q'(r,v,t)  , 


and 


=  Kfn-i(r,v,t)  . 


By  hypothesis 


(19) 


(20a) 


(20b) 


(21) 


Furthermore 


\|ri(r,v,t)  =  KQ'(r,v,t) 

^  Mo’maxW  ^  •  (22) 

Continuing  the  iteration  we  find 

'|fn(r,v,t)  s  M(c"„^anax)"  •  (25) 

Thus,  the  Neumaimseries  converges  pointwise,  and  the  theorem  is  proved.  The  fact 
that  \|((r,v,t)  is  positive  follows  from  the  fact  that  every  term  in  the  series  is 
positive. 
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Next,  assume  that  c"(r,v)  Is  not  bounded.  Then  we  can  state: 

Theorem  2.  Let  Q'(r,v,t)  be  bounded.  Then  a  unique,  positive  solution  to 
the  transport  eqtiatlon  exists  which  is  a  continuous  function  of  t  and  r 
and  an  Integrable  function  of  v. 

Proof:  The  proof  proceeds  along  the  lines  of  Theorem  1;  however,  it  is  suf¬ 
ficient  to  show  that  the  series 


^^3v\|fn(r,v,t) 


(2i+) 


converges  polntwlse.  Consider  first 


^^vfj  (r,v,t)  S  mJ  dt'J^  )  ,v'-*v)  .  (25) 

However,  we  have  assumed  in  Sec.  Ill  that  a{r,v'-»v)  could  be  written  in  the  form 


a(r,v'-^)  =  ^Ni(r)ai(v'-*v)  , 


(26) 


where  the  Ni^(r)  are  bounded,  say  by  Niq.  Then 


y  d^vt,(r,v,t)  S  Mt  v'd^'6i(v' )ai(-' 


(27) 


But  li(v')  and  ai(v')  are  assumed  in  Sec.  Ill  to  be  bounded,  say  by  and 

a i o .  Then 


^d®vi|(,(r,v,t)  i  v'd®v'  , 


(28) 


and  since  v'  S  (cf.  Sec.  Ill),  we  have 
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d®VT|/j^(r,v,t)  S  Kt  ,  0  S  K  < 


(29) 


Similarly,  we  find 


d  v\|/j^(r,v,t)  S 


n  I 


,  0  «  K'  < 


(50) 


and  the  theorem  is  proved.  Again  each  term  is  positive,  so  that  the  solution  is 
positive. 

Next  consider  the  case  that  Q'(r,v,t)  is  not  bounded,  but  integrable.  The 
following  theorems  are  simple  to  prove  by  straightforward  construction  of  the 
Neiunam series,  as  above. 

Theorem  3*  If  Q'(r,v,t)  =  Qo(£»Z)®('*')>  where  Qo(£>v)  is  bounded,  then, 

for  the  time-dependent  transport  equation 

(a)  If  c"(r,v)  is  bounded,  a  unique,  positive  solution,  \|((r,v,t), 
exists.  \|/(r,v,t)  is  a  continuous  function  of  r  and  v,  and 
itr(r,v,t)-Q' (r,v,t)  is  a  continuous  fimction  of  t. 

(b)  If  c"(r,v)  is  not  bounded,  then  the  solution,  i|((r,v,t),  may  be 
an  integrable  rather  than  a  continuous  function  of  v.  Otherwise 
the  conclusions  are  unchanged. 

Theorem  4.  If  Q'(r,v,t)  =  Q,(r,t)6(v-Vo),  where  Q,(r,t)  is  bounded,  then, 

(a)  If  va(r,v'-*v)  is  boxinded,  a  unique,  positive  solution  \|r(r,v,t) 
exists.  \|f(r,v,t)  is  a  continuous  function  r  and  t,  and  t(£>v,t)- 
Q'(r,v,t)  is  a  continuous  function  of  v. 

li^ 


(b)  If  va(r,v'-^)  Is  not  bounded,  then  i|r(r,v,t)-5[(r,v,t)  may  be  an 
Integrable,  rather  than  a  continuous,  function  of  v*  (in  proving 
part  (b)  we  make  use  of  the  fact  (cf.  Sec.  Ill)  that  v  and  have 

both  been  assumed  bounded.) 

Theorem  5*  If  =  Qj2(v,t)8(r-rQ), where  Q'2(r,t)  Is  bounded,  then 

(a)  If  c"(r,v)  Is  bounded,  a  positive,  unique  solution  \|f(r,v,t) 
exists.  \|f(r,v,t)  Is  a  continuous  function  of  v  and  t,  and  an  Integr¬ 
able  function  of  r. 

(b)  If  c"(r,v)  is  not  bounded,  then  \|;(r,v,t)  may  be  an  Integrable 
rather  than  a  continuous  function  of  v. 

The  theorems  for  the  cases  In  which  Q'(r,v,t)  involves  a  delta  function  in 
more  than  a  single  variable  may  easily  be  constructed  by  appropriately  combining 
the  above  theorems.  We  shall  not  state  them  separately. 

V.  The  Time -Dependent  "Adjoint"  Equation 
Here  we  deal  with  Eq.  (lO). 

The  theorems  will  all  be  stated  without  proof,  since  the  proofs  are  completely 
analogous  to  those  given  in  the  previous  section. 

Theorem  6.  If  Q'(r,v,t)  Is  bounded,  then  a  continuous,  \inlque,  positive 
solution  exists.  (Note  that  In  proving  this  theorem  it  is  necessary  tc 
make  use  of  the  fact  that  c(r,v)  Is  bovmded  (cf.  Sec.  III).) 
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Theorem  7.  If  Q'(r,v,t)  =  ^(r,v)8(t),  where  ^  is  bovinded,  then  a 

unique,  positive  solution,  \|r(r,v,t)  exists.  \ii'(r,v,t)  is  a  continuous 
function  of  r  and  v,  and  >if(r,v,t)-5'(r,v,t)  is  a  continuous  function 
of  t. 

Theorem  8.  If  Q'(r,v,t)  =  Qi(r,t)6{v-Vo),  where  5^  is  bounded,  then 

(a)  If  va(r,v'-*v)  is  bounded,  a  vinique  positive  solution,  \|^(r,v,t) 
exists.  \if(r,v,t)  is  a  continuous  function  of  r  and  t,  and  \ir(r,v,t)- 
5E'(r,v,t)  is  a  continuous  function  of  v. 

(b)  If  va(r,v'-*v)  is  unboimded but  |^'(r , v)  bounded,  then  the  con¬ 

ditions  above  hold  except  that  \i?(r,v,t)-Q'(r,v,t)  may  be  an  integrable, 
rather  than  a  continuous,  function  of  v. 

Theorem  9-  If  §'(r,v,t)  =  Qa(v,t)8(r-ro),  where  §2 is  bounded,  then  a 

a  unique,  positive  solution,  ii?(r,v,t),  exists.  \ii'(r,v,t)  is  a  continuous 
function  of  v  and  t,  and  an  integrable  function  of  r. 

Again  when  5'  Involves  delta  functions  in  more  than  a  single  variable,  the 
appropriate  theorems  can  be  constructed  by  appropriately  combining  the  results 
above,  and  so  we  shall  avoid  stating  them  explicitly. 

VI.  The  Time-* Independent  Transport  and  Adjoint  Equations 

In  certain  rather  restrictive  cases,  it  can  be  proved  that  unique,  positive 
solutions  of  the  time-independent  transport  and  ajolnt  equations  exist.  In  gen¬ 
eral,  the  restrictions  are  much  more  severe  than  is  the  case  for  the  time-dependent 
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equations;  ve  shall  see  later  that  this  Is  to  be  expected. 


Theorem  10.  Let  Q'(r,v)  be  botinded  and  positive.  Then  If  c'(r,v)  <1,  a 
unique,  positive,  and  continuous  solution  to  the  time -Independent  trans¬ 
port  equation  (Eq.  (?))  exists. 

We  note  that  the  Integral  equation  (7)  differs  from  Eq.  (lb)  In  that  In  the 
former  the  dependent  variable  Is  the  angular  flux,  ♦(rjv),  rather  than  the  angular 
density  \^(r, v) .  Thus  all  theorums  which  we  shall  prove  Involving  the  time- 
independent  transport  equation  may  not  apply  to  the  angular  density  unless  \|f(r,v) 
vanishes  sufficiently  rapidly  as  v^.  This  Is  not  a  real  worry,  since  one  Is 
generally  Interested  In  the  flux  rather  than  the  angular  density. 

As  usual,  we  prove  the  theorem  by  constructing  the  Neumann  series : 

<'(r,v)  =  ,  (31a) 

with 


^oir,y)  «  Q'(r,v)  , 


(31b) 


and 


<'n(£>Z)  =  A“<i(r,v)  *  A''’n-l(l>l)  •  (51c) 

By  hypothesis,  ♦o(£>Z)  bounded  and  positive: 

0  i  ♦o(£'Z)  «  M  <  »  .  (32) 
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Next  assume  bounded  by  M',  say.  Then 


♦n(l»l)  *  A<'n-j^(r,v) 

S  M’c'niaj^  dRexp 


|^a(r,r-Rn,v^ 


a(r-Rfi,v) 


(33) 


The  integral  is  easily  shown  to  be  equal  to  unity.  Thus 


♦n(r,v)  S  M'c'max  <M'  ,  (3»t) 

since  we  have  assumed  c'(r,v)  <  1.  This  proves  the  theorem  since  the  Neumannseries 
converges  pointwise,  and  each  term  is  seen  to  be  positive. 

The  next  theorem  is  readily  proved  in  essentially  the  same  manner. 

Theorem  11.  If  Q'(r,v)  =  Qj^(r,v)6(v-VQ),  where  Qj^(r)  is  bounded,  then 

if  c'(r,v)  <1  and  a(r,Vo-^)  is  bovinded,  a  unique,  positive  solution,  \|f(r,v), 
exists.  t(r,v)  is  a  continuous  function  of  r  and  \|f(r,v)-Q' (r,v)  is  a 
continuous  function  of  v. 


Theorem  12.  If  a(r,v)Q' (r,v)  is  Integrable ,  then  if  c(r,v)  <  1,  a  unique, 
positive,  integrable  solution  exists. 

This  is  the  theorem  proved  by  Olhoeft  (2) .  The  procedure  is  to  construct  the 
Neumam series  for  the  collision  density  x(£>£)  defined  by 

x(r,v)  =  a(r,v)*(r,v)  ,  (35) 


and  prove  that 
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y'd®vd3rl^n(£>l)  ^  J  d®rd®vXn(r,v) 


(35) 


Since  the  details  are  given  in  Ref.  2,  they  will  be  omitted  here^ 

We  next  turn  our  attention  to  the  time-independent  adjoint  equation.  Here 
the  theorems  are  quite  analogous  to  theorems  10-12  proved  for  the  transport 
equation,  except  that  the  roles  of  c  and  c*  are  reversed.  Since  the  proofs  are 
so  similar,  we  simply  state  the  theorems. 

Theorem  13.  Let  i  5!'(r,v)  be  bounded  and  positive.  Then  if  c(r,v)  <  1, 

V - 

a  unique,  positive,  continuous  solution  exists  to  the  time -independent 
adjoint  equation. 

Theorem  l4.  Let  i  =  Qi(r)8(v-VQ) ,  where  Qi(r)  is  bounded  and 

positive.  Then  if  c(r,v)  <  1  and  o(£,X*Zo)  is  bounded,  a  unique,  positive 

solution,  ?(r,v),  exists.  t(r,v)  is  a  continuous  function  of  r  euid 
~  1  ~ 

ilr  -  —  Q  is  a  continuous  function  of  v. 

Theorem  15.  If  c'(r,v)  <  1,  and  if  i  a(r,v)§'(r,v)  is  Integrable,  then  a 

~  V 

unique,  positive,  integrable  solution  exists. 

Note  that  in  each  of  these  theorems  we  have  placed  restrictions  on  —  Q' 

V 

rather  thsin  on  Q' .  If  we  make  the  reasonable  assumption  that  no  zero-energy 

source  neutrons  are  introduced  into  the  system,  then  the  conditions  can  equally 

well  be  stated  as  conditions  on  Q' rather  than  on  i  Q'. 

V 

Next  we  prove  a  uniqueness  theorem  for  cases  in  which  existence  has  not  been 
proved. 
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Theorem  l6.  If  a  solution  of  class  L  of  the  time-independent  adjoint  equa¬ 
tion  exists,  then  aui  solution  to  the  time -independent  transport  equation, 
if  it  exists,  will  be  unique. 

We  sketch  the  proof;  Suppose  there  are  two  solutions  to  the  transport  equation, 
\1(^  and  \|fg.  Then  \|f  =  obeys  the  eqxxation 

v(n«7+c(r,v))il((r,v)  =  J  v'a(r,v'-*2)'l'(r,v' )d^'/  (37) 

with 


'l'(rg,v)  =  0  ,  n’no  <  0. 


(38) 


Consider  the  adjoint  equation  with  zero  outgoing  angular  density. 

v(-n*7-t-a(r,v)  )\if(r,v)  =  q(r,v)-^vo(r,V'+v' )ii'(r,v' )d\' ,  (39) 


^(rsiv)  =  0  ,  n'no  >0  • 


(i^o) 


(We  have  hypothesized  that  such  a  solution  exists.)  If  we  now  multiply  Eq.  (37)  by 
t(r,v),  and  multiply  Eq.  (39)  by  '|f(r,v),  subtract  and  integrate  over  r  and  v  we 
obtain 


J 


daly/  d^'nQi|/(r,v)\if(r,v) 


d^'v'a(r,v'-»y) 


X  i|f(r,v')ii((r,v) 


-  J d®v'va(»,V'+v' )+(r,v' )t(r,v) 

-  f  <i^v\|f(r,v)q(r,v)  .  (4l) 


20 


The  left  hand  side  of  this  equation  has  been  obtained  vlth  the  help  of  Gauss's 
Theorem, and  by  vlrture  of  Eqs.  (58)  andOt^O),  It  vanishes.  Similarly,  the  first 
two  terms  on  the  right  side  of  Eq.  (4l)  cancel,  and  we  obtain 


y'd^r^ 


d^v\|^(l,y)q(x,y)  =  0 


Since  q(r,v)  Is  arbitrary  (and  positive)  It  follows  that 


(42) 


+(i:,y)  s  0  , 


(43) 


proving  the  theorem. 

We  note  that  the  conditions  of  this  theorem  are  satisfied  for  (among  other 
cases)  c(r,v)  <  1;  thus,  Eq.  (39)  can  posess  no  solution  for  c(r,v)  <  1  subject 
to  the  boundary  conditions  (38).  Since  these  are  precisely  the  equations  satisfied 
by  the  neutron  density  In  a  reactor,  we  have  succeeded  In  proving,  as  a  by-product, 
the  not  surpylslng  result  that  a  reactor  cannot  be  critical  If  fewer  neutrons  are 
emitted  than  absorbed  following  each  collision  (c  <  l).  A  simpler  proof  of  this 
theorem  (and  of  the  following  theorem)  for  the  case  that  only  scattering  and  pure 
absorption  are  present,  (l.e.,  no  fission)  Is  given  In  the  Appendix. 

The  analog  of  Theorem  (16)  Is; 

Theorem  IJ.  If  an  L  solution  of  the  time-independent  transport  equation 

exlsti^  then  an  L®  solution  of  the  time -independent  adjoint  equation.  If 

It  exists,  will  be  unique.. 

The  proof  of  this  theorem  Is  essentially  identical  with  that  of  Theorem  16. 

The  Implication  of  this  theorem  Is  that  a  reactor  cannot  be  critical  If  c'(r,v)  <  1. 
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This  can  be  seen  from  a  slight  modification  of  the  arguments  used  above  to  shov  the 


same  result  for  c(r,v)  <  1. 


VII.  Some  Miscellaneous  Results 


In  the  previous  section  we  have  pointed  out  that  the  uniqueness  of  the  time- 
independent  solution  is  equivalent  to  the  statement  that  a  reactor  cannot  be  critical. 
Thus  we  have  shown  (Theorem  l6)  that  no  reactor  can  be  critical  for  o(r,v)  <  1, 
which  is,  of  course,  physically  obvious.  However,  such  a  condition  Is  somewhat  too 
stringent  since  it  is  clear  that  for  a  system  of  finite  size  neutron  leakage  may 
prevent  a  reactor  from  becoming  critical  even  for  c  somewhat  larger  than  one. 

We  can  obtain  an  estimate  of  the  minimum  value  of  c  for  which  a  reactor  of  a 
given  size  will  become  critical  from  Theorem  12.  The  proof,  which  we  omitted  In 
Sec.  VI,  Involves  the  construction  of  the  Neumann  series  solution  to  the  following 
Integral  equation  for  the  collision  density  x  *  o*: 


x(r,v)  =  a(r,v)Q'{r,v)+Ax(r,v)  , 


(1^4a) 


where  A  Is  the  integral  operator 


Af( 


r,v) 


dfVv'-^d^' 


r'exp 


■a(r,r',v)|  o(r,v)8(r ' -r+RJl)  ^  f(r',v')  . 

__  __  _  _  _  -  - 


(44b) 


In  the  proof,  one  shows  that  If 


JJ  d®vd®rXji(r,v)  S  M  , 


(45a) 
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that 


J  J  d3vd3rxn+j(r,v)  S  Cm,vM  •  ('»5b) 

Actvially,  a  somewhat  stronger  condition  holds,  since  from  Eq.  it  follows  im¬ 

mediately  that 

vl  )tf(r,v'H^)-  1  , 

-*  a(r,v) 

m 

In  obtaining  the  exponential  in  Eq.  (U6)  was  set  equal  to  zero.  If  instead 

we  approximate  Eq.  (^^6)  by  the  following  expression 

x^{r,v')i^Td^'  ,  (4?) 

we  see  that  an  approximate  limit  for  the  convergence  of  the  Neuman  series  is 

c  <  .  (48) 

1-e-^ 

Here  t  is  the  average  optical  chord  length  of  the  system  and  c  is  the  average  value 
of  c(r,v) . 

This  expression  should  give  a  rough  estimate  of  critical  size  since  one  can 

write 


J J d®vd®rAXn-n(l<z)  =  d^'(l“exp  "0(i:>r+Rsn> 


T  m  'St  , 


(49) 


where  a  is  the  average  cross-section  and  t  is 


the  average  chord  length (“47/8). 
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The  averages  of  c  and  a  are  taken  both  with  respect  to  r  and  v. 

A  second  resxilt  ^Ich  follows  from  the  theorems  proved  In  the  earlier  section 
Is  that  the  solution  of  the  time -dependent  transport  or  adjoint  eqiiatlons  can  In¬ 
crease  no  faster  than  exponentially  If  the  source  Is  of  exponential  order  or  less. 
Consider  first  the  transport  equation.  Let  us  write  the  source  In  the  form 


Q"(r,v,t)  =  '  Qe(r,v)f(t)  . 


(50) 


For  simplicity  we  shall  assume  that  and  c"  are  bounded.  Then  we  wish  to  examine 
the  time-dependence  of  the  Neumann  series,  Eq.  (19).  Since  Q  Is  bounded  we  have 


i  Mf(t)  , 

S  Mc"niaxOmax^  tf(t)dt  , 


(51) 

(52) 


and  In  general 


dt  jX  •  •  •  ^  )dtn 


’•'n  *  ”(®"max°max 


(55) 


Now  from  Euler's  Identify 


/'k-r’v(tn)dt„  =  r^IlLVr!.(V)dt. 

i  'i  (n-i); 


(5M 


we  have 


+n  *  **(®"max'^max^  ^max  ~ 


(55) 


Thus 


2k 


(56) 


♦  S  Mfaax 


00 


n»o  n.' 


or 


t  *  ^max®*P(‘^"n!ax°max^^  • 


(57) 


Thus  the  time  behavior  of  \|f  is  asymptotically  given  either  by  f(t)  or  by  the 
c"  0  "t 

exponential  e  ,  whichever  is  more  important  at  large  times. 

If  f(t)  a  e  ,  the  integral  in  (5^)  can  be  evaluat6d  explicitly. 


In(t) 


e“°*^dT 


(58) 


Clearly 


In(t) 


(n-ijr  3at  a 


)  , 


(59) 


so  that 


,  -Oft 


*n  S  K- 


(^"maxc^max)  (-1)“  l-e 


(n-i)!  9a’ 


,n-i  (  Q  )  > 


(6o) 


for  n  >0.  Then 


\|r  s  M  {l+c"inaxffmax®'^  Y  i£j5S*^S2Ll*^  (-l)"^  —  (i^^  ))/^^^^ 

^  n.'  9oP  a 

n=o 


or 
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max°niax 


where  Taylor's  theorem  has  been  used  to  perform  the  sum.  We  note  that  this  limit 
is  always  positive. 

A  very  similar  argument  can  be  used  to  delimit  the  time  behavior  of  the 
adjoint  equation;  we  shall  omit  it  here. 


Appendix 


Theorems  (16)  and  (l?)  nay  be  proved  in  a  different  manner  from  that  used  in 
the  text  If  we  assxame  that  no  mechanism  is  present  for  neutron  regeneration  other 
than  ordinary  elastic  and/or  inelastic  scattering.  (We;.  note  that  in  this  case 
o{r,v)  SI.)  In  such  a  case,  we  expect  that  the  kernel  a(r,v’-^y)  obeys  the 
principle  of  detailed  balsmce 

v'M(v' )a(r,y '■*y)  =  vM(v)a(r,y+y' )  ,  (A-l) 

since  In  the  absence  of  sources  and  sinks  of  neutrons,  an  equilibrium  should  be 
approached.  The  equilibrium  spectrum,  M(v),  Is  the  Maxwell -Boltzmann  distribution 

I 

Consider  now  Eq.  (37)  of  the  text,  with  the  boundary  condition,  Eq.  (38)* 
Define  a  new  dependent  variable  by  the  relation 

♦(r,y)  =  F(r,y)\|f(r,y)  ,  (A-3) 

where  P(r,y)  Is  some  well-behaved  but  otherwise  arbitrary  function  to  be  chosen 
later.  Then  Eq.  (37)  becomes 

V.  V  1  +  =  r v'o(v'-^)  »  (A-^) 

- F  F  J  -  -  F(y') 

with 

♦(5>v)  -  0  ,  O'no  <  0  .  (A-5) 


M(v)  ~ 


2 

V  exp 


-87 


(The  function  F  must  be  chosen  so  that  £q.  (A-^)  Is  still  satisfied.) 

If  Eq.  is  multiplied  by  ♦(r,v)  and  Integrated  over  r  and  v  we  obtain, 

with  the  aid  of  Gauss's  Theorem 


♦^(r,v) 

F(r,v) 


+ 


d  wa(r,v) 


F(r,v) 


(A-6) 


Next  consider  the  identity 


[♦(rjv)  -  'l'(r,v') )  £  0  , 


(A-7) 


or 


2  2 
♦  (r,v)'>(r,v' )  S  ^  (r,v)  +  <'(r,v')  } 


(A-8) 


Then,  by  vlrture  of  (A-9,  the  right  side  of  Eq.  {A-6)  may  be  written 


f  f  d^vd®v'v'a(r,v'-^) 

JJ  --  -  F(r,y') 


i  f  /d^d^v'v'  '*S)  {>^(r,v)+*^(r,v' ) )  ( 

2J  J  F(r,y') 

i  [  d®v'v'a(x,x')c(x,x')  [  dVv'V 

F(r,y')  F(r,v') 


(A-]D) 
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If  we  now  choose  F(r,v') 


1 

M(v') 


the  second  term  on  the  right  of  Eq.  (A-IO) 


becomes 


d®vd®v'v'M(v')cr(r,v'-nr)»^(r,v) 


=  ^ J  d®vd^'vM(v)a(r,v^' )*^(r,v) 


(A-11) 


by  Eq.  (A-l).  This  becomes,  upon  Integration  over  v', 


I  ^^wM(v)o(r,Y)c(l,Y)ti^(x,Y) 


(A-12) 


writing  M(v)  for  l/r(x,v)  throughout  Eq.  (A-6),  and  using  (A-9)>  (A-10),  and  (A-12), 
we  find 


/ds/  d®vnQ*v<'^(r,v)M(v)  +  fi^rf  d®'vv^(v)a(r,v)'l>^(r,v) 


d®wa(r,v)c(r,v)vM(v)<t>^(j;,v)  S  0  . 


(A-13) 


The  first  term  on  the  left  side  of  Eq.  (AJ.5)  is,  by  virture  of  (A-5  ) ,  non-negative. 
The  other  two  terms  combine  to  give 


d^vvM(v)a(r,v)<'^(r,v)[l-c  (r,v)  1 


(A-14) 


which  is  always  non-negative  since  c(r,v)  <1.  However,  Eq.  (A-13)  tells  us  that 
the  sum  of  these  terms  is  non-positive.  Thus 


4'(r,v)  s  0  , 


(A-15) 


29 


and  the  theorem  is  proved.  We  see  that  a  reactor  with  no  fuel  can  never  be  critical. 


The  uniqueness  of  the  time-independent  adjoint  equation  for  the  care  of  no 
reproduction  can  be  proved  in  an  entirely  analogous  fashion.  This  method  is  at 
slight  generalization  of  that  used  in  Ref.  (l)  for  proving  uniqueness  of  the  one- 
speed  equations. 
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